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Abstract. - We consider a quantum system of fixed size consisting of a regular chain of 
n-level subsystems, where n is finite. Forming groups of A'' subsystems each, we show that the 
\ strength of interaction between the groups scales with A*'~ As a consequence, if the total 

. system is in a thermal state with inverse temperature /?, a sufficient condition for subgroups of 

size to be approximately in a thermal state with the same temperature is y/N ^ P6E, where 
SE is the width of the occupied level spectrum of the total system. These scaling properties 
indicate on what scale local temperatures may be meaningfully defined as intensive variables. 
This question is particularly relevant for non-equilibrium scenarios such as heat conduction etc. 
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- ^ , Introduction. - Thermodynamics has successfully been applied to the description of 

' macroscopic systems since more than a century. For that reason this theory enjoys widespread 
^ I acceptance. Nevertheless its microscopic foundation is, in most cases, not well understood. 

The theory of thermodynamics is based on the notions of extensive and conjugate intensive 
thermodynamic variables. The existence of the thermodynamical limit in a mathematical 
sense has been shown for many fundamental cases 012113 • The standard proofs are based 
on the idea that, as the spatial extension increases, the surface of a region in space grows 
slower than its volume. If the interaction potential is short-ranged enough, one can show that 
the "effective" interactions between one region and another become negligible in the limit of 
infinite size, implying extensivity. 

However, the scaling behavior of the interactions with the size of those regions has, to our 
knowledge, not been studied in any detail yet. While it has been shown that they vanish in 
the thermodynamical limit, it remains unclear what happens in situations that are in some 
sense only close to this limit. There is not even a precise understanding of what "close" means 
in that case. 
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For standard applications of thermodynamics this might not pose any serious problem since 
the number of particles within any region is so large that deviations from infinite systems may 
safely be neglected. Nevertheless the differences should become important when the considered 
regions are significantly smaller. It is here where a quantum approach becomes imperative 

The applicability of thcrmodynamical concepts to mesoscopic or even microscopic systems 
has intensively been discussed in recent years [Hj, since nano-scale physics has attracted in- 
creased attention j^j. Due to the advance of experimental techniques the measurement of 
thermodynamic quantities like temperature with a spatial resolution on the nanometer scale 
seems within reach These techniques have already been applied for a new type of 

scanning microscopy, using a temperature sensor jlOIll) . that shows resolutions below lOOnm. 
An important question thus arises 0: Down to what spatial scale does a meaningful notion 
of temperature exist at all? 

In this paper we consider a quantum system of fixed Hilbert-space dimension dim, com- 
posed of identical elementary n-level subsystems. We form identical groups of N subsys- 
tems each {dim — n"'^^) and show that the interaction strength between the groups scales 
inversely proportional to \/N. Based on this scaling relation and assuming the total system 
to be in a thermal state we analyse, for what group size N a thcrmodynamical description 
of the individual group is appropriate. The results are confirmed by numerical studies for a 
chain of 8 spins. For this chain partitions into two 4-spin-groups, four 2-spin-groups and eight 
single spins are considered. 

Scaling Law. - We consider a chain of identical n-level systems with identical nearest 
neighbour interactions. The Hamiltonian of such a linear chain may be written as 

H = J2 '^'°^(^') + '^"^tijj + 1), (1) 

3 

where the index j labels the elementary subsystems. The local terms hioc{j) and the nearest 
neighbour interactions hint{j,j + 1) have the form |12j : 

hiocU) = ^ X! ^"(-^^ 

a 

2 

h^nt{j,J + 1) = ^ I] <Ja{j) (j + 1) (2) 

a, 13 

Here, the aa{j) are SU{n) generators with a, /3 = 1, 2, . . . , — 1, n being the dimension of 
one subsystem. For the SU (n) generators we adopt the trace relations ■ 

Tr[a„W] = 

dtTTl 

TT[aa{i) (Tpij)] = 2——SijSap (3) 

where the trace (Tr) has to be taken over the whole system of dimension dim. As a consequence 
the trace over one elementary subsystem reads: Tr^- [cTa{j) cr^O')] = 2 (5^/3- The coefficients in 
equation (O are then given by Aq, = Tt{H a^ii)) / dimandCap ~ Tt{H a^ii) (Tjsii + l)) /dim, 
respectively, and taken to be independent of i. We assume periodic boundary conditions and 
H to be traceless. 

If we now form uq groups of N subsystems each, we can split the Hamiltonian H into two 
parts, 

H = H^j^^ + In (4) 
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where H^"^ is the sum of the Hamihonians of the mdividual groups and In describes the 
interaction between each group and its nearest neighbour, 



We label the eigenstates of the total Hamiltonian H and its energies by greek indices and 
eigenstatcs and energies of the group Hamiltonian -ff^^ by latin indices, i. e. 



H\fi)=E^\fi) and H^°^ Ij^) 



(6) 



Obviously, the jj^v) are simply products of group eigenstates. We now proceed to compare 
two characteristic quantities of the system according to equation Q). 
First, consider the identity 



\{JnW^ ^ {]N\I%\jN) 



(7) 



where the left hand side can be interpreted as the mean squared energy distance between a 
level E^ and the energy Ej^^. The average of equation ^ over all is identical with I^, 



where 



dim 



(8) 



a, 13 



Here we have used equation 0. 

The second quantity of interest is the width of the distribution of the energy levels of the 
total system 5E around the mean energy E ~ TT{H)/dim = 0: 



6E 



1 1 2 


Tr 


\h-e)'' 


dim \ 


dim 



I — n 



Q,/3 Q 



(9) 



Here we have again used equation (|3J)- Combining equations ||HJ) and (Q, we get the following 
scaling law: 

\n<^JE (10) 



TV 

Note that this law is a property of the Hamiltonian |^ and does not depend on the state 
of the system. The equality sign applies, if the local terms vanish, = 0. Equation 1)10(1 
does not depend on the convention used in equation since normalisations cancel. It is 
straightforward to see that the same scaling law holds for a three dimensional lattice with 
cubes of subsystems as the subgroups. 

Relation H10|) has been derived here for a perfectly homogenous system. However, it 
is evident by means of stationary perturbation theory that it still holds approximately for 
systems with small disorder. 

Distributions. ~ We now use the scaling law ()l()|l to estimate the density matrix elements 
of the groups of N subsystems assuming that the total system is in a thermal state with the 
density matrix 



it^lpW) 



(11) 
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in the eigenbasis of H . Here, Z is the partition sum and (3 — {ksT) ^ the inverse temperature. 
Transforming the density matrix into the eigenbasis of h'"^^ we obtain 

(Mj) = E lO'l-")!' = E IOV>P cxp(-/3(£;^-i?,)) (12) 

/J A' 

(For simphcity we skip the index N from now on). 

A thermodynamical description of the groups of N subsystems by canonical density ma- 
trices with the same inverse temperature (3 is appropriate, if was approximately pro- 
portional to exp(—/3 iSj). Truncating the sum to terms with f3\E^ — £'j| <C 1 (which will be 
motivated below), we can expand the second exponential in the rhs of equation (|12|l into a 
Taylor series around Ej up to second order, 

m3) « ^ - mm + ^(ji/'ij) j , (13) 

where we have used J2fj,^ft = + (jI^Ij) and equation 0. The second order term 

has to be taken into account, since the sum over all energies E'^ adds up positive and negative 
contributions for the first order terms but only positive ones for the second order terms. Since 
(jI-^I.?)^ (il^^li): both correction terms are small if 

P^Mm«l■ (14) 

Equation H14|l justifies the above truncation if the distribution Wj{ii) ~ times the 

density of states ry(-E'^j — Ej) decayed faster than exp(— /3|£'p — E-j\) for \E^ ~ E-j \ > y/ 
and fixed j. We numerically verified this behavior for a class of systems as shown in figure 
El If, on the other hand, the truncation of the sum was not possible, the rhs of equation 1)1 3|l 
would contain additional correction terms invalidating a local thermodynamical description. 

In the basis the off-diagonal elements of the density matrix exp(— /?i?^) 

vanish for \Ej — Ej>\ > 21 because \{j\^J'){^^\j')\ ~ 0. When \Ej — Eji\ < 21, one can use the 
same approximation as for the diagonal terms, where now the zero order term is zero. In 
the first and second order corrections each term of the sum carries a phase and thus these 
corrections are smaller than for diagonal elements. 

Combining equations (|10|l and H14|) we thus conclude that the condition 

VN:^(3 5E. (15) 
is sufficient to allow for an approximate local thermodynamical description for a group size N. 

_____ _2 

In equation H15|) we have used that I is the arithmetic mean of all these are positive 

quantities and therefore /3I ^ 1 implies (3y/ <C 1 for almost all states Equation 
i|15|l is the main result of our paper. 

Numerical results. - To test the condition ^5(1 and the approximations involved, we 
investigate a chain of 8 spins with a Hamiltonian of the form |^ rewritten as 

AE * 8 3 

= — E '^-(j) + ^ E E ^"'^ ^a(j) apij + 1) (16) 

J = l 3 = 1 0,0=1 

Without loss of generality, we restrict the local part to terms in only; different local terms 
would merely imply a rotation of the coordinate system. Periodic boundary conditions are 
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chosen, <7a{^) = o'ai^)- A is a scaling factor for the interactions and the ctq (a — 1,2,3) are 
the Pauli matrizes. In the following all energies (i.e. A, (3^^) will be taken to be in terms of 
AE. 

For each realisation the matrix elements c^^ are randomly chosen from the interval [—1, 1] 
with equal weight. Based on 100 such realisations, we find the average of 5E to be « 5A. 
Therefore, from equation (|15|) . a local thermodynamical description is expected to be appro- 
priate for VN I {PX) > 5. 

We consider three different partitions of this system into groups of adjoining subsystems as 
described in equation Q). The partitions are: two 4-spin-chains, four 2-spin-chains and eight 
single spins. For the inverse temperature of the total system (3 and the interaction strength A 
we consider the values PX = 0.1, 0.2, 0.3, 0.4 which implies 2.5 < ^/N / (/3A) < 20. 

First we test the scaling behavior of the interaction, equation (|10|) . In figure we have 
plotted the ratio l/SE (sec equations (0 and©) for the three partitions and for 
comparison. The result confirms equation (|10|l . The values arc slightly smaller than l/ViV, 
which is due to the neglect of the local terms of the Hamiltonian. 

Then we calculate Wj{p,) rj{Ei^i — Ej) for all three partitions. For each j we have plotted the 
distribution Wj{iJ,) rj{E^~Ej) versus x = {E ^ — Ej) / ^J{j\P\]) thus rescaling its width to unity. 
Figure 121 shows a superposition of all theese plots for one realisation. For comparison we have 
plotted in the range, where \E^ — Ej\ > ^ (jV^Ij), the functions 0.25 exp(— — Ej\) = 
0.25 exp(-0.5|a;|), taking |3^/{J\P\j) = 0.5 < 1. The normalisation 0.25 was deliberately 
chosen to show that the distribution Wj{p,) ri{Ef^ — Ej) indeed decays fast enough in that 
range. All numerical tests we made showed such a behavior. 




Fig. 1 Fig. 2 



Fig. 1 - The ratio of the average interaction strength and the level spreading 1/5E (dots) as a function 
of the number of subsystems per group A'^ for A = AE. The line shows 1/^/N. 

Fig. 2 - The conditional probabilities \{fJ-\j)\'^ weighted with the density of states ri[E^) as a function 
of X = - Ej)/^{j\P\j). The lines 0.25 exp(-0.5 are shown for comparison. 

Let us now assume the whole system to be in a thermal state with a reciprocal temperature 
(3 and a density matrix according to equation (|11|) . We calculate the distance between this 
density matrix and a product of canonical density matrices of the subgroups corresponding 
to a partition 

p^W^prij) (17) 

where PQ^'ij) is of the same form as in equation but only for one subgroup. As a measure 
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Fig. 3 - Distance between the density matrix p and a product of canonical density matrices of the 
subgroups corresponding to the selected partition /5. A'^ is the number of subsystems per group. 

Fig. 4 - Ratios of the spectral inverse temperatures Pspec of the subgroups to the given inverse 
temperature /3 of the total system. A'^ is the number of subsystems per group. 



for the distance we use 



dtstip,p) = ^Tv[ip-pr] (18) 

The result is shown in figure |31 The distance between the state of a global and a local 
thermodynamical description is found to be approximately proportional to /3 A and decreases 
as the group size N increases. The points with \/N / (/3A) > 5, for which our estimates should 
apply, are below 0.1, which supports equation (|15|l . 

Finally, to further confirm our findings, we calculate a "spectral temperature" of each 
subgroup in the following way |18) : The canonical density matrix of the whole system is 
transformed into the product basis ©. To each pair of states formed by an excited energy 
level Ei and the ground level Eq of a subgroup, one can then assign a Boltzmann factor with 
an inverse temperature /?'. The spectral temperature is the sum of all these f3' weighted by 
the occupation probability pi of the excited level Ei: 

Pspec- 2^i_ E E ^^^^ 

Such a /3spec can be defined for any state and coincides with the thermodynamical (3 for a 
canonical state. Since periodic boundary conditions were assumed, all subgroups of the same 
size N have the same temperature. In figure 0] the ratios between fispec defined by equation 
l)19|l and the inverse temperature of the total system (3 are plotted versus group size N . Here 
the points, which fulfill our criteria, lie above 0.9, which again confirms our results. The 
spectral temperatures of the subgroups are, in general, higher than the temperature of the 
total system. This indicates that the state of a subgroup is more mixed than the state of 
the whole system, which may be due to the entanglement of the subgroups caused by their 
mutual interaction |14) . 

Conclusion. - We have studied a modular system of fixed size composed of elementary 
subsytems with a finite energy spectrum and nearest neighbour interactions. We have shown 
that if one forms groups of N subsystems each, the interaction between neighbouring groups 
scales as N~ We have then considered a chain of such interacting subsystems in a global 
thermal state with canonical density matrix. We have argued that due to the scaling property 
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of the interaction, the reduced density matrix of each group may be approximated by a 
canonical one with the same temperature as that of the total density matrix, if ^/N ^ (3 SE: 
The temperature becomes an intensive quantity on a coarse-grained size and length-scale only. 
In the same way, energy becomes more and more extensive as the group-size increases. 

We have tested this assertion numerically with a chain of 8 interacting spins and find that 
our predictions are met in spite of the still small size of the total system. 

Our studies should be extended to the total system being in a local thermodynamical 
equilibrium only. Here, heat conduction [BJ^] becomes an interesting problem. For strong 
coupling; meaningful temperature-profiles can be defined with limited resolution only 

* * * 

We thank M. Michel, H. Schmidt, M. StoUsteimer, and F. Tonner for fruitful discussions. 
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